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We introduce two families of sum-of-squares (SOS) decompositions for the Bell operators associ¬ 
ated with the tilted CHSH expressions introduced in Phys. Rev. Lett. 108, 100402 (2012). These 
SOS decompositions provide tight upper bounds on the maximal quantum value of these Bell ex¬ 
pressions. Moreover, they establish algebraic relations that are necessarily satisfied by any quantum 
state and observables yielding the optimal quantum value. These algebraic relations are then used 
to show that the tilted CHSH expressions provide robust self-tests for any partially entangled two- 
qubit state. This application to self-testing follows closely the approach of Phys. Rev. A 87, 050102 
(2013), where we identify and correct two non-trivial flaws. 


I. INTRODUCTION 

We consider a bipartite Bell scenario [1] with two bi¬ 
nary observables Aq, Ai for Alice and two binary ob¬ 
servables Bq, Bi for Bob, such that =1, By = I, and 
[Ax, By] = 0 for all x, 2 / = 0,1 [2]. The following family of 
tilted CHSH operators for this scenario was introduced 
in 0: 

la = ChAq -\- AqBq + AqBi -\- AiBq — AiBi , (1) 

where a G [0, 2[ is a parameter and a = 0 corresponds 
to the CHSH operator. One of our objectives in this 
article is to derive an upper bound rja on the maximal 
quantum expectation value of Iq,, that is to show that 
(Tq) < Va for all possible quantum states and measure¬ 
ment operators Ax and By . This is equivalent to showing 
that the operator = rjai — la is positive semidefinite, 
i.e. la cz 0, for any measurement operators Ax and By. 
This in turn can be proven by providing a set of opera¬ 
tors {f)} which are polynomial functions of Ax and By 
such that 

Ta=Y.P^P^ (2) 

i 

holds for any set of measurement operators satisfying the 
algebraic properties = I, = I, and [Ax, By] — 0. 
Indeed, if Xa is of this form, it is obviously positive 
semidefinite. A decomposition of the form ([^ is called a 
sum of squares and can be defined in full generality for 
other Bell operators than the ones considered here. 

Though it may be difficult to find a SOS decomposition 
for an arbitrary Bell operator, especially one yielding a 
tight bound r]a on the maximal quantum violation, once 
one has been found, verifying that ([^ holds, and thus 
that (la) < Vcn usually involves only a few simple com¬ 
putations. That is, a SOS provides a simple certificate 
that (Xa) < rja- Furthermore, the search for optimal SOS 
can be cast as a series of semidefinite programs (SDP) 
that turns out to be simply the dual formulation |4] of 
the SDP hierarchy introduced in BM- Finally, as shown 
in [7], an optimal SOS, i.e. one for which (la) < rja is 


a tight bound, provides useful information about the op¬ 
timal quantum strategy and can find an application in 
robust self-testing. 

Self-testing is the process through which one can guar¬ 
antee that two devices satisfy certain properties, e.g. that 
they implement measurements on a quantum state which 
is close, up to a local isometry, to a given reference state, 
only by observing the correlations in a Bell experiment 
[5]. The possibility of self-testing means that the inter¬ 
action on a classical level with quantum devices can be 
sufficient to assure users that they indeed hold devices 
that conform to an ideal specification: as long as the 
self-testing criteria are satisfied, it is guaranteed that the 
devices have not malfunctioned or been tampered with. 
Effectively, this allows the users to treat their devices as 
black boxes. This is the core idea of device-independent 
quantum information processing, where self-testing has 
been used as a primitive to establish schemes for verified 
quantum computing and cryptographic tasks such as 
randomness expansion mm- 

Despite their interest, very few examples of explicit 
SOS decompositions for Bell operators have been given 
in the literature. One example is a SOS for the fam¬ 
ily of “guess your neighbour’s input” inequalities intro¬ 
duced in m- A second example is a SOS for the CHSH 
operator appearing in [^, which we recover in this arti¬ 
cle. Finally, SOS decompositions for the entire family of 
tilted CHSH operators ([^ were proposed in [7]. How¬ 
ever, as we point out here, these SOS only hold in the 
range a G [0.156,1.955]; outside that range they are not 
valid SOS due to a sign error that cannot be simply fixed. 
Furthermore, even within the validity range, these SOS 
decompositions are not sufhcient for the self-testing ap¬ 
plication proposed in [7] (see Appendix [C|) . 

In this article, we introduce two different simple SOS 
decompositions for the tilted CHSH operators ([^ that 
are valid for the entire range a G [0,2[, including the 
CHSH case a = 0. Our systematic approach to this 
problem is of independent interest as it can probably 
be adapted to other Bell operators. We then show, fol¬ 
lowing [12] and [7|, how to apply these SOS decompo¬ 
sitions to robust self-testing of any partially entangled 
state. Moreover, we take the self-testing analysis further 
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by showing that the isometry providing a robust self-test 
for the state, also provides a self-test for the action of 
the measurement operators. Incidentally, we point out 
and fix a small mistake in m concerning the regulariza¬ 
tion procedure used to define the local isometry used for 
self-testing. 

This article is organized as follows. We first present 
our SOS decompositions and detail the approach we used 
to obtain them. We then sketch their application to 
self-testing, which is comprehensively discussed in Ap¬ 
pendix]^ The flaws that we identified in [7] and m and 
how our work resolves them is presented in Appendix [C| 


II. SOS DECOMPOSITIONS FOR TILTED 
CHSH INEQUALITIES 

We first start by reviewing some of the properties of 
the tilted CHSH expressions We then introduce the 
concept of SOS decompositions and derive such decom¬ 
positions for the tilted CHSH operators. 


A. Optimal quantum strategies for 

The optimal quantum value for the Bell expressions 
0 was computed in [3] by optimizing explicitly over all 
quantum states and measurements. The optimal value is 
jmax _ ^8 _|_ 20 ^, to be compared to the classical bound 
2 -I- a. An interesting property of this class of Bell op¬ 
erators is that for all values of a G [0, 2[, the maximal 
quantum value can be achieved by a partially entangled 
qubit pair 

I')/’) = cos 0 |00) -I- sin0 |11) (3) 

with a = a{6) = 2/vT+~2~taj?~^, effectively covering 
the full range of partial entanglement in two-qubit states, 
9 G ]0,7r/4]. The operators used to achieve this maximal 
quantum violation are 

Aq = cr^ Ai = (4a) 

Bq = cos sin B\ = cos ^Oz — sin ^ ax (4b) 

where tan fi = sin 29 and ax^z are the x, z Pauli operators. 


B. SOS decompositions for Bell operators 

Our goal is to find a decomposition as in ([^ in terms of 
a set of polynomials {Pi} when the constant term in Iq, is 
the maximal quantum value of the tilted CHSH operator, 
Va = For simplicity, we restrict our search space 

to the span of a small canonical basis of nine monomials 
of degree two, 

(5) 


To simplify the search, we will later pick a different basis 
of this vector space of polynomials, {Ri}i. Thus, the op¬ 
erators Pi are decomposed as Pi = (using implicit 

summation on repeated indices). The expression of la 
as an SOS ([^ then becomes 

ia = ^ • (6) 

If q} is now seen as the fth component of a vector the 
hermitian matrix M is then the Gram matrix of this set 
of nine vectors and is therefore positive semidefinite. The 
converse is also true: if M is positive semidefinite, there 
exists a (non-unique) set of vectors {q'^} such that the 
components are the hermitian products (q'^jq'^) = 
gf and therefore any operator of the form of the RHS 
of 0 with M ^ 0 is a sum of squares. A set of such 
vectors {q^'} is given by the columns of any matrix square 
root of M, such as its Cholesky decomposition if it is 
nonsingular. 

We are now looking for a positive semidefinite matrix 
M such that 0 holds. This equation imposes linear 
constraints on M. Indeed, we can decompose both sides 
of the equality Iq = in a basis of the quadratic 

products of all elements in which is of size 25 

rather than 81 due to the algebraic relations satisfied by 
the measurement operators. A canonical basis for these 
products is 

^i+AB = {I; ^1; ^i^o} 

'S'{i, Bq, Bi, BqBi, BiBq} ■ (7) 

Writing = FJ^ibRi where Ei runs through 

and each A* is a matrix of complex coefficients, and like¬ 
wise la = s^Ei, the SOS condition 0 reduces to 

s* = Tr(M^J^*) (i = I,...,25). (8) 

We are thus left with a set of 25 linear equality con¬ 
straints on M as well as the positive semidefiniteness 
constraint AI A 0. This is a semidefinite programming 
feasibility problem, which can be approached with nu¬ 
merical tools such as SeDuMi m- 

However, while we could attempt to recover the exact 
analytical expression that the numerical solution approx¬ 
imates, this will not be a good approach for a continuous 
class of Bell operators. We will thus tackle this problem 
analytically, and this requires that we simplify the prob¬ 
lem as much as possible. One first simplification comes 
from our knowledge of a state and measurements which 
achieve the quantum bound. One such system is specified 
in Eqs. 0 and 0, consisting in the partially entangled 
qubit pair \ijj) = cos 9 |00) -|-sin0 |I1) with spin measure¬ 
ments along given axes. Because this strategy achieves 
the quantum bound the expectation value of 

vanishes. As a consequence, any SOS decomposition for 
Xa of the form 0 must have each of its terms vanish in 
expectation as well. Hence, a valid SOS decomposition 
for la must be made up of terms for which Pi\ip) = 0 
in this maximally violating quantum system. This last 


Sl+AB — {I, Aq, Ai} ® {I, Bq, Bi} . 
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equation defines four constraints that all Pi must obey, 
one per basis vector of the Hilbert space. Indeed, writing 
the most general P in our search space as r • V where 

V = (I, Hq, Hi, Hq, BijAqBq, AqBi, AiBq, AiBi) (9) 


and demanding that the four components of P li/i) = 
P{cos6 |00) +sin0 |11)) (with the observables specified in 
Q) vanish, we find four independent linear constraints 
on the vector r. We thus find that the space of candidates 
Pi is spanned by the following five operators: 


Zx — Zi'Q , (10a) 

I - , (10b) 

cXa — sZaHTb — Wa^b , (10c) 

cXb - sXaZb - Za^b , (lOd) 

sXa^b — ZaZb + cZa , (lOe) 


where c = cos 20, s = sin 20, and the Z and X operators 
are defined as 


Za — Aq 


Bq + Bi 
2cos^ 


Xa — Ai , 


Bq — Bi 
2 sin /i 


(lla) 

(llb) 


It is easily verified that these five operators indeed vanish 
in the explicit two-qubit system we considered, as Z and 
X are then the Pauli operators. 

This first step lets us write the Pi polynomials as linear 
combinations of five operators instead of nine, hence the 
M matrix in our SDP feasibility problem is now 5x5. 

One further simplification comes from exploiting the 
symmetry of the tilted CHSH operator to impose a sim¬ 
ilar symmetry on its SOS decompositions. The effect of 
symmetries on SOS decompositions has been studied in 
the case of polynomials of commutative variables in El, 
and the following takes inspiration from those results. We 
observe that changing the sign of Ai -A —Ai and swap¬ 
ping Bq o Bi leaves la (and therefore la) invariant. 
This transformation induces a representation of the cyclic 
group C 2 on the vector space of operators span(iSi+AB). 


Let a G End(span(5i+AB)) be the endomorphism rep¬ 
resenting the transformation. Due to the invariance of 
the tilted CHSH operator, the five-dimensional subspace 
spanned by operators (10) is itself invariant under a by 


definition. Therefore, we separate this subspace as a di¬ 
rect sum of irreducible representation subspaces, all of 
dimension one according to the representation theory 
of cyclic groups. We then choose our basis operators 
{Ri}i=i to be the basis elements associated with such a 
decomposition. This means that a acts on this basis as 
(j{Ri) = ±i?i, depending on i. 

We now apply this symmetry transformation to both 
sides of ([^: 


Ta = cT{Ia) = = rIm'^'^R^ , ( 12 ) 


where M'^'' = depending on the sign brought by 

the transformation of the basis operators. This means 
that M' represents a new SOS decomposition for Xa- 
Moreover, because convex combinations of SOS decom¬ 
positions are also valid SOS decompositions of the same 
operator, a third SOS decomposition is found: = 

(M -|- M')/2. This last SOS has the property of being 
invariant under the symmetry transformation. 

While asymmetric SOS may still exist, this result is 
useful because it lets us focus on symmetric decompo¬ 
sitions, which have a smaller number of degrees of free¬ 
dom. Indeed, because elements in M' and M only differ 
by their sign, the symmetrization of M to will take 
to zero the elements that change sign. The elements in 
question correspond to indices (p,, v) such that and 
Ri, span representation subspaces of different irreducible 
representations of € 2 - As a result, the symmetrized SOS 
matrix is block-diagonal, with one block associated 
to each of the two irreducible representations. 

Considering the discussion thus far, we now choose a 
basis for the subspace containing the SOS polynomials. 
The basis we choose is {i?i = • V} where the vec¬ 

tors are defined as such (we label the columns with the 
operators defining V for convenience): 


ri 

r2 

rs 

U 

rs 


I 

( » 
( 

^ Vl + 


Aq 

-2 

\/l + s'^ 
0 

0 

0 



0 


Ai Bq Bi AqBo AqBi AiBq AiBi 

0 1 1 0 0 0 0 ), 


0 0 0 1 


1 0 



0 

-2 

Vl + 

-2c 

Vl + 


0 0 

0 0 

1 -1 


c c 

1 -1 

0 0 


1 

c 

1 



(13) 
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It is easily checked that these basis operators separate the 
space in two isotypical subspaces, i.e., subspaces that fall 
under the same irreducible representation of the cyclic 
group: -Ri, 2,3 are invariant under the symmetry transfor¬ 
mation of la, while i? 4_5 change sign. The block structure 
of symmetric SOS matrices is therefore 3 0 2, where the 
first block corresponds to the trivial representation and 
the second to the parity representation where the group 
generator is represented by — 1 : we have cr(i?i) = Ri for 
i = 1, 2, 3 and a{Ri) = —Ri for i = 4, 5. 

We now examine the problem of finding SOS decom¬ 
positions with a = 0, where Xq is the CHSH operator, 
and with a = Q!( 7 r/ 8 ). 



FIG. 1. Solution set for symmetric SOS decompositions for 
the CHSH operator in Si+ab- 


1. The CHSH case: a = Q 


We hrst look at the determination of SOS decompo¬ 
sitions in the simplest case, where a = 0. Without yet 
taking symmetry into account, the linear constraints im¬ 
posed on M by Eq. Q imply (omitting the lower trian¬ 
gular part given by hermiticity) 


M = 


2V2 


/A 0 0 5 0 \ 

j q 0 —5/2 
t 0 -5/2 

fj, 0 

q . 


(14) 


with (7 = 1 — A — 7 and t = A 0 7 — /r. Thus, M de¬ 
pends on four parameters, three of which are real due to 
hermiticity. The fourth parameter, 5, can be taken to 0 
by the symmetry argument given above. We see that, in 
this specihc case, the 1 0 2 0 1 © 1 block structure of a 
symmetric M goes further than what is imposed by the 
symmetry of ■ This is not surprising, as the CHSH op¬ 
erator, lacking a marginal term, is more symmetric than 
the general la operator. This makes it straightforward 
to impose positive semidefiniteness: a necessary and suf¬ 
ficient condition is the positive semidefiniteness of the 
blocks. We find the conditions 


A,/i,9,t>0, (15) 

Xt-q^>0. (16) 


We therefore end up with a simple description of the full 
solution set of symmetric SOS decompositions of Xq in 
up to the computation of a square root for M to 
get the actual polynomials Pi in (|^. We represent this 
solution set in terms of A, /i and q in Figure 


In this solution set, the extremal points are the most 
interesting as they generate the entire set by convex com¬ 
binations. Here, the extremal points are the five vertices 
of the set and the smooth quadratic surface resulting 
from inequality (16). Because of their situation at the 


intersection of inequality constraints, the M matrices at 
the hve vertices are of low rank, which eases the deter¬ 
mination of a square root. Those distinguished extremal 
points are listed in Table |Ij We note that the first four 
vertex SOS in the table are equivalent up to multiplica¬ 
tion on the left by a dichotomic operator on one or two of 
the operators that are squared in the sum. For instance, 
the identity R 2 = AqRi implies that R\ = which 
means that the SOS Ci and C 2 are equivalent. 

We reproduce here two SOS decompositions of Xq re¬ 
sulting from those extremal points: 


Xq = 


1 


4v^ L' 


Io + 2(ZaXb0Xa^b)' 


(17) 


Io = ^ [(-^a 0 ^b)" + (-Xa 0 Xb)"] (18) 


where = ^ 0 , = ^ 1 , Zb = {Bq + Bi)/\/2, and 

Xb = {Bq — Bi)/'j2. We note that the decomposition 
(fl^, which we denote as C 4 in Table III also appears in 

r 

Two additional SOS decompositions for CHSH are also 
given in Appendix [B| 


2. The 6 — 7r/8 case 

Next, we choose the (arbitrary) value of 0 = 7 r/ 8 . The 
equality constraints (|^ give to the SOS matrix the fol¬ 
lowing form: 




M = 


8^2 


-1 

V2 


3(5 



P + h 

i -27 

-y/25 




2 0 37 - A 

^5 

-25 

(19) 




V 
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FIG. 2. Solution set for symmetric SOS decompositions for 
2 ^a(7r/8) in Si+AB- This set is delimited by a quadratic surface 
(invariant along /3) and a cubic surface. Yang and Navascues’ 
SOS [7] is represented as a red dot. 


where, again, (5 can be set to zero by symmetry to give 
M a block-diagonal structure. 

Contrary to the CHSH matrix, the block structure in 
this case is coarser with fewer, larger blocks, which results 
in more nonlinear positive semidefiniteness constraints. 
The positive semidefiniteness of M is checked by using 
a generalized form of Sylvester’s criterion, namely that 
a matrix is positive semidefinite if and only if all of its 
principal minors (submatrices obtained by jointly elimi¬ 
nating lines and columns of the same indices), including 
the matrix itself, have a nonnegative determinant m- 
The solution set is visualized as a three-dimensional con¬ 
vex region in Figure delimited by a quadratic surface 
and a cubic surface. 

We identify on Figure two points that stand out, 
namely the two cusps at the intersection of the two sur¬ 
faces delimiting the region. Their coordinates are found 
to be (/3,7,A) = (1/2,-1/3,0) and (3/2,1/3,8/3). Both 
blocks of M have rank one at these points, which eases 
the computation of square roots N that describe the SOS 
polynomials. We find respectively 


N2 


Ni 


73 / l /72 -1 

872 1 0 0 


-1 0 0 \ 

0 0 -Ij ’ 



( 20 ) 

( 21 ) 


The rows in these matrices give us the coefficient vectors 
Qi for two explicit SOS decompositions ofIa{-K/8)- 


3. General tilted CHSH 

The general approach is the same as above: first, the 
linear constraints ([^ are imposed, and imposing sym¬ 
metry leaves us with three real degrees of freedom. The 
positive semidefiniteness of the blocks of M is then en¬ 
forced by applying the generalized Sylvester’s criterion 
m- This amounts to one cubic, four quadratic and five 


linear inequalities. For a given value of a, the solution 
set is readily visualized as a convex set in 

In the previous section, we examined the case of a = 
a(7r/8) and identified two extremal points in the solution 
set where each block of the SOS matrix M is of rank 
one. When we carry out the explicit expansion of the 
sums of squares, these decompositions hint at an easy 
generalization to all values of a for both points. This 
leads us to the two SOS decompositions that follow. Let 
us first define the following permutations of the CHSH 
operator: 


5' = Ho(Ho-Hi)+Hi(Ho + Hi), (22) 

5" = Ho(Ho+Hi)-Hi(Ho-Hi), (23) 

S'” = ^o(^o — Bi) — Ai{Bq + Bi). (24) 


The two decompositions we find are the following: 

1 


T = 

rv 


2B 


xl + {aA^-S'f 


(25) 


and 


T = 

rv 


2B 


{2Ao - C‘ 

+ {2Ai — 


.Bn 


+ Bi a „ 
2 ^ 2^^ 


Bq — Bi 


a 


S” 


(26) 


The two SOS decompositions given in the CHSH case 
0 , ^ are limiting cases of J2^ an d (|2^ . Indepen¬ 
dently of the way we arrived at (25) and (26), it is readily 
verified that they correspond to valid SOS. 


III. APPLICATION TO SELF-TESTING 


We conclude by discussing the application of the above 
results to self-testing. We know that measuring the state 
I')/) = cos9 |00) -I-sin6* |11) with the observables @ leads 
to the maximal value of the Bell expression la ■ We 
say that this Bell expression provides a robust self-test 
for this particular reference state and reference observ¬ 
ables if the converse also holds, in a noise-tolerant way — 
that is, if an expected value of la close to the maximum 
jmax necessarily corresponds to measurements involving 
a state and observables that are close to the reference, 
up to a local isometry. More precisely, we show, fol¬ 
lowing the framework of [12] and [7|, the following. Let 
{la) be the expectation value of the Bell expression la 
obtained by measuring a physical state lip) with physi¬ 


cal observables A^ and By. Let 


and Ax and By be 


the reference state and reference observables correspond¬ 
ing to th e optimal quantum strategies for defined in 
Then if (1^) > — e, there exists a local 


HA 


Sec. 

isometry $ = ® '1 >b and a state |junk) such that 


||4>(A2, (g) By Ipj)) - iJunk) (g) {Ax (g By) \ip)\\ < e' (27) 

for x,y € {—1,0,1}, where the subscript —1 refers to the 
identity operator and where e' = 0{y/e). The precise 
relation between e and e! is given in Appendix [A| 
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FIG. 3. The local isometry used to locate the partially en¬ 
tangled qubit pair in the physical system. The dashed line 
shows the physical separation between Alice (top) and Bob 
(bottom). 


The case of the CHSH inequality (a = 0) was already 
considered in [T^] , where Eq. (27) was obtained but with 
a weaker e' = robustness. Equation in the 

case a;, y = —1 for arbitrary values of a € [0, 2[, corre¬ 
sponding to a self-test of any partially entangled state 
(but not of the associated observables) was considered in 
[7], where the authors showed that the use of SOS decom¬ 
positions could in principle lead to an optimal e' = 0(-y/e) 
robustness. This conclusion, however, does not follow 
from the analysis presented in j7] due to different short¬ 
comings in the derivation of intermediate results (see Ap¬ 
pendix]^ for details). Here we resolve these shortcom¬ 
ings, thanks in particular to the two SOS decompositions 
(25) and (26), and establish the robust self-testing con¬ 
ditions thus furth6r Bxtcnciiiig the cuicilysis in j7] to 
the self-testing of the observables (|^ in addition to the 
state ([^. 

Following closely [7], we now sketch the proof of the 
robust self-testing result (27) and refer to Appendix [A| 
for details. The proof proceeds in three steps. First, 
we define as in |12j the isometry appearing in (27) as the 


successive action of a set of gates, represented as a circuit 
in Figure]^ acting on the initial state |00) 0 |^). The 
gates are defined in terms of the measurement operators 
Ax and By as an inversion of ([^ in an attempt to recover 
the behaviour of the Pauli operators. Specifically, we 
define as in (11) 


Zp^ — Aq 

Bq + Bi 


Zb = 


Xa = A,, 

Bq — Bi 


2 cosy 


Ab = 


2 sin y 


(28a) 

(28b) 


Because we need unitary gates in the isometry, a regular¬ 
ization procedure is applied to normalize the eigenvalues 
of the operators: we define Z'^ = Z^\Z^\~^ where Z^ 
is Zb with its zero eigenvalues changed to 1. The same 
procedure is applied to define Ag. Alice’s operators are 
unchanged because they are already unitary: Z'a = Za 
and = Xa- As a side-note, we point out a slight 
mistake in |12| concerning the above regularization pro¬ 
cedure, see Appendix [C| 

Second, we observe, as mentioned earlier and as 
pointed out in [7], that the existence of a SOS decom¬ 
position for la of the form implies that any state \ij)) 
and operators A^ and By achieving the quantum bound 


jmax obey the relations A |V^) = 0. Moreover, this 
property is robust: if instead the expectation value is 
(Iq) = — e, then we have ||A |'0)|| < -s/e- This ob¬ 

servation applied to the two SOS (25) and (26) implies 
the following identities up to robust error terms: 


{Za-Zb)\^)=0, (29a) 

(sin6»AA(I + Zb) - cos6»Ab(I - Za)) IV’) = 0, (29b) 

{ZaXa - XaZa) I^) = 0. (29c) 


For instance, identity ( 29a[ ) in the case a = 0 is easily 
seen to follow from the term Pi = {—Za + Zb) present 
in the SOS decomposition (18). The three identities 


for arbitrary a follow in a similar, though more involved, 
manner from the two SOS (25) and (26), see Appendix[A| 
The third step of the proof then examines the action 
of the isometry on a state and observables satisfying the 
algebraic relations (29). It establishes that in this case 
the relations (27) necessarily hold between the inputs and 
outputs of the isometry. As in [7], a lengthy series of 
triangle inequalities is needed to go from the SOS bounds 
on ||7i |V’)|| to the robustness bounds for (29) and finally 
to (27). We refer to Appendix|^for the derivation, where 
explicit bounds in 0(-\/e) for all the identities in (27) are 
determined. 


IV. DISCUSSION 

In this work, we provided tools to simplify the search 
for sums-of-squares decompositions for Bell operators, 
exploiting their symmetries and the knowledge of systems 
that maximally violate the associated Bell inequalities. 
We applied this approach to find two SOS decomposi¬ 
tions for the family of tilted CHSH inequalities of [3] (as 
well as some extra decompositions in the special case of 
the CHSH inequality). We made use of these new SOS 
decompositions to complete and extend the proof of Yang 
and Navascues [7] by showing that a close-to-maximal 
quantum violation of a tilted CHSH inequality provides 
a robust self-test for the reference state and measurement 
operators associated with the inequality. 

The general form C-^/e of our distance bounds is opti¬ 
mal in the noise parameter e in the sense that a larger 
exponent for e would contradict the hypothesis (Iq,) < e. 
These distance bounds, though, become very sensitive to 
noise as the entanglement diminishes. This is not sur¬ 
prising, as taking the entanglement parameter 9 close to 
zero takes the extremal quantum behaviour for the Iq, 
inequality closer to the local set. 

The reference systems for the tilted CHSH inequalities 
are particularly relevant to randomness generation, as ex¬ 
plained in [3]. Our self-testing statement may therefore 
be useful in establishing device-independent protocols us¬ 
ing partially entangled states as a resource for random¬ 
ness expansion. 

We hope that a similar approach to finding SOS de¬ 
compositions for Bell inequalities will find further appli- 
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cations in robust property testing in different systems, 
for instance with higher-dimensional reference Hilbert 
spaces. 
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Appendix A: Robust self-test 


1. Introduction 

In this section, we show how our SOS results can be 
used to provide a robust self-test of the systems that 
maximally violate the Xq, inequalities. We work in the 
following setting: two isolated parties, Alice and Bob, 
hold black boxes that each take one of two inputs, re¬ 
spectively X and y, in {0,1}, and each return one of two 
outputs, a and 6, in {-fl, —1}. They wish to make sure 
that their boxes share the partially entangled qubit pair 
|'0) = COS0 100) -|-sin0 111) and that the observables char¬ 
acterizing their measurements on their share of the state 
are given in terms of the Pauli operators as 

Ao = 0-2 Ai= (Ta; (Ala) 

Bq = cos y az + sin yox Bi = cos y az — sin y Ox (Alb) 

with tan/r = sin 2d. In reality, the boxes hold the phys¬ 
ical state I'f/j) and the observables are the hermitian and 
dichotomic (i.e. of eigenvalues ±1) operators Ax and By 
- this is general because there is no assumption made on 
the dimensionality of the Hilbert space, which can thus 
be extended to purify the state and make the measure¬ 
ments projective. 

We follow the framework of McKague et al. m and say 
that a Bell expression provides a robust self-test for a ref¬ 
erence system if, for a violation of the corresponding Bell 
inequality that is e-close to the quantum maximum (we 
say that the system satisfies the self-testing criterion), 
the state and measurements that give rise to this viola¬ 
tion are close to the reference, with the error vanishing as 
e goes to zero. Formally, this is stated as the existence of 
a local isometry with respect to Alice and Bob that takes 
!'(/>) and the action of their observables Ax and By to a 
system close to the reference ji/i) and By, in tensor 
product with uncorrelated degrees of freedom. The pre¬ 
cise mathematical statement is the following: assuming 
that the self-testing criterion holds, there exists a local 
isometry $ = G and a state |junk) such that 

\\<^{Ax®By !■!/())-|junk)(g)(A,c(g)Hy) \'ip)\\ (e) (A2) 


for all X, y S {—1,0,1} (where operators with subscript 
— 1 refer to the identity), and lime_>o exy(e) = 0. This 
is understood as meaning that there exists a procedure 
(that need not be accessible experimentally) that Alice 
and Bob can follow locally to perform a change of basis 
after which the reference state and operators (or some¬ 
thing close) are found in a four-dimensional subspace of 
the global Hilbert space. 

We show in this section a robust self-test for the 
ideal reference system highlighted above for Xq, for any 
a G [0, 2[, with self-testing criterion {tp\Ia |V') = — e 

(with Xq, expressed in terms of Ax and By), leading to 
explicit self-testing bounds Cxy G 0{^/e). This result uses 
the same techniq ues a s [7] to achieve the self-test of the 
state (statement ( |A2[ ) for x = y = — 1) and extends the 
reasoning to the self-test of operators. 

Following [3 [12], we defin e the following operators, 
essentially an inversion of (Al I for the physical operators: 


Za — Aq 

Bq + Bi 
2 cosy 


= 


= Al, 

Bq — Bi 
2 sin y 


Ar = 


(A3a) 

(A3b) 


We regularize these operators into unitaries to use them 
in the isometry. The procedure is as follows: start from 
the physical observable, for example and change all 
zero eigenvalues to one, resulting in a new hermitian op¬ 
erator Z*. Then, normalize the eigenvalues by defining 
Zg = ZgjZgl”^. This last operator is by construction 
unitary, self-adjoint, commutes with Zb, and has the 
property that Z^Zb = I-^bI- Similarly, we define Ag, 
Z^ and X'p^, noting that Alice’s regularized X'p^ and Z^ 
operators actually coincide with the physical operators 
Xa and Ab. 

We now define the self-testing isometry by representing 
it in circuit form in Figure In general, the action of 
the isometry on the physical state |'0) is 


+ Ai(I-Zi)(I + Z{i)|^)|10) 

+ A{i(i + z;)(i-z{i) 1^)101) 

+ AiAB(I-^;)(I-^B)l^)|ll)]- 


(A4) 


In this form, the motivation for this isometry is readily 
understood. In the reference system with state \tjj) and 


operators (Al), the A and Z operators above (and their 
regularizations) are the Pauli operators ax and az ■ Then, 
$(|^/i)) is easily seen to be |00) \tj}), where the partially en¬ 
tangled qubit pair |i/i) has been extracted onto the ancilla 
qubits, leaving behind the state |00) in the physical reg¬ 
ister, unentangled with the ancilla. Moreover, by manip¬ 
ulating (A4) using the operators’ properties, it is easily 


seen in this ideal setting that acting on the input with 
the measurement operators is equivalent to acting on the 
corresponding ancillae with the reference operators. 

Our goal is to extend this statement to all other sys¬ 
tems that maximally violate the Xq, inequality, and show 







its robustness to noise by deriving (A2). To do this, we 


need to show that the self-testing criterion implies some 
robust properties on the action of the physical operators 
on the state. 


S 5 = i?3 (and S5 = ra), which is then expressed in terms 
of the other four Si as 


—c 


S 5 =Ai\ —S 2 - 


\/r 


^4 . 


(All) 


2. Application of our SOS decompositions to 
robust self-testing 

We first start by reformulating our two SOS decompo¬ 
sitions and show that, even though they contain only four 
independent terms, they still provide a way to test the 
entire five-dimensional space of relations that we defined 
as we identified candidates for SOS decompositions in 
Eqs. (10). Defining the following polynomials in Ax,By: 


Si=I<, 

52 = aAi - S' 

5 3 = 2 Ao - ir 

5 4 = 2Ai - /r 


, Bq -|- Bi a 


S" 


. Bq — Bi a 


2 ^ 


with 


S' = Aq{Bq — Bi) + Ai{Bq + Bi ), 
S" = Aq{Bq Bi) — Ai{Bq — Bi) , 
S'" = Aq{Bq — Bi) — Ai[Bq + Bi) , 

the two SOS decompositions of are then 


(A5a) 

(A5b) 

(A5c) 

(A5d) 


(A6) 

(A7) 

(AS) 


T — 
— 


1 


2 /" 


:(^f + 5|) = 


1 


2 /" 


:{Sl + Sl). (A9) 


We can also express the Si polynomials as • V with 
V defined in ([9) and expressed in terms of the basis 
vectors (13): 


51 = cri - r2 - ra 

52 = -rs , 

sa = -VT+^ri -I- 


S4 = 


-2 


vT" 


: r4 


vTT7 

c 


:(r2-ra). 


vT~ 


rs. 


(AlOa) 

(AlOb) 

(AlOc) 

(AlOd) 


The four Si polynomials above do not span the en¬ 
tire candidate subspace span({i?i}) that we identified in 
Eqs. (10). Indeed, while S2 and S4 generate the parity 


isotypical subspace spanned by r4 and rs, the other iso¬ 
typical subspace is not spanned by the other two vectors. 
For example, ra cannot be decomposed in terms of Si and 
S3 alone. We note however that left multiplication by Ai 
takes an operator from the parity to the identity repre¬ 
sentation subspaces and vice versa. In fact, the simple 
relation R 3 = A 1 R 4 holds. Because R 3 is not a linear 
combination of the Si operators while R 4 is, we define 


Suppose now that the physical state \ip) and observ¬ 
ables Ax, By satisfy the self-testing criterion with error e, 
i.e., the expectation value (Zq) for the Xa Bell expression 
is such that (2^^) = — e. Denote as Xa the Bell oper¬ 

ator defined in terms of the observables A^ and By, and 
let {Pi} be a set of polynomials in these observables such 
that /““I —Xa=Xa= X)i Pi Pi- Then, the self-testing 
criterion directly implies the bounds {-ipl P- < e 

or, equivalently, ||Pi |'0)|| < ^/£■ Hence, given a sum- 
of-squares decomposition of Xa, the action of the observ¬ 
ables on the state is constrained for any state satisfying 
the self-testing criterion. 

The SOS decompositions found in this article are valid 
for any A^ and By that are hermitian and dichotomic, 
where A^ and By commute. By their definition and the 
assumption that Alice and Bob are separate, the physical 
observables A^ and By match these properties. Hence, 
we can derive from our two SOS decompositions some 
useful bounds on the action of the observables on the 
state. 


In our SOS decompositions (A9), the Pi operators are 
(2/™®'X)“i/^S'i. We thus let S = y^2l2^-\/e so that the 
self-testing criterion implies \\Si \ip)\\ < (5 for z = 1, 2,3,4. 
From this, we deduce a similar bound that we will need 
for the action of S5. Because Ai is unitary and appears 
as a left multiplier in S 3 , we can drop it from the norm 
IIS's I'0)11 to find the norm of a linear combination of S2 
and S4 acting on |0). We then use the triangle inequality 
to bound this by 0 (i/e): 


11^5 10)11 = 




< 211^2 10)11 


■01 -|- ■ 


\/e 


-S4I 


1^4 10)11 


C -I- VT+02 
< -;- 0 . 


(A12) 

(A13) 

(A14) 


Now that we have 0(5) = O(0e) bounds on \\Si |0)|| 
for all i, by using the triangle inequality we can robustly 
certify the action of the entire subspace of operators in 
Si+AB that vanish in the ideal setting, which proves to 
be very useful in showing the effectiveness and robust¬ 
ness of the self-testing isometry. Furthermore, we can 
extend this space beyond Si+ab by left-multiplication 
by bounded operators as was done for S3, which we will 
use in the self-test of measurements. 

Recalling our general SOS decompositions, we notice 
that the four Si polynomials in addition to the S3 poly¬ 
nomial (all formulated in terms of A^ and By) can be 
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linearly combined to form the following operators, used 
in [7] to prove the self-test: 

Za - Zb = . (A15) 

sind^A -I- — cosdXe 

= (<' “ ^ (AM) 

In anticipation to the measurement self-test, we add to 
this list the anticommutator between Alice’s two observ¬ 
ables. This polynomial is not in 5 i+ab, but we can use 
left-multiplication to decompose it in terms of Sf. 


ZaXa + XaZa — 


Vl + s^ 


2 S 2 -|- Aq 


+ Ai 


(^ — cS2 + \/l -f S^S' 4 ^ 

(^cSi + \/1 


(A17) 


We note that the anticommutator {.^bj A'b} is zero from 
definition (A31. We will also need the following, which 
we write in terms of for the sake of briefness: 


Ab - -1a(I - cZb) = . (A18) 

s 2 s 

We now write down the explicit robustness bounds de¬ 
rived from the relations above in the same fashion as 

(lAni): 


show, the regularization procedure applied to define the 
unitaries in the isometry only introduces new error terms 
of the same order of 0 (^/e) as the bounds for the unreg¬ 
ularized operators. 


a. Self-testing the state 


We focus in this section on the claim of [7], that is, the 
self-testing bound (A2) for the state (i.e., x = y = —1). 

We start by showing an analogue of inequality (A19) 
for the regularized operators. To do so, we use the trian¬ 
gle inequality to separate \\{ZIq — Z'a) \ip)\\ into two terms. 
The first one is bounded as follows: 


IK^b-^b) 1^)11 


= \\zi,(i-\ZB\)m 

(A27) 

= 11(1-IZa^bI) 1^)11 

(A28) 

< 11(1-ZaZb) 1^)11 

(A29) 

VI 

(A30) 


In the first equality, we used the identity Z^\Z^\ = z^ 
that we highlighted when defining the regularized oper¬ 
ators. In the second equality, we use the unitarity of Z'q 
and the fact that the absolute value of an operator is 
unchanged by acting on the left with a unitary opera¬ 
tor, here Za- The l ast in equality uses the unitarity of 
Za again to recover (A19). Recalling that = Za, we 
thus have 


||(Za-^b) 1^)11 <<5i, 
^sin QXa (I -I- ^b ) — cos 0 Ab (I — ^A)) 


(Ab-s-IAa(I-cZb)) 


< <^5 I 


\\{ZAXA + XAZA)m\<5-, 


A 
a 5 


with 


r /I \ C T V^I -(- ^ 

= (I + c) . (5 


<52 = 

<55 = 




2s2 

((2-c) + \/r + ^)5, 


Vl + 

8sind 
Vl + (c -I- y/l -I- s^) 


4s 


yr 


2s2 


-I- c -I- 1 / 1 - 1 - 52 ^ i5 . 


(A19) 

< <52 ) 
(A20) 

(A2I) 

(A22) 

(A23) 

(A24) 

(A25) 

(A26) 


We now prove that the isometry defined earlier pro¬ 
vides a robust self-test for our reference system. This 
consists in applying a series of transformations on bounds 
(A19)-(A22| in order to reach bounds (A2). Concretely, 


in order to prove the self-testing bounds, the action of 
the unitary operators that constitute the isometry must 
be shown to be restricted by the self-testing cri terion . 
If there were no need for regularization, bounds (A19|- 


(A22) would directly apply to the isometry. As we will 


ll(Z'-Zi) 1^)11 

<\\{z'^ -z^)m + \\{z^ -z’^)m (A31) 

< 2<5i. (A32) 


We then prove a bound for Ag similar to ( A30| ) , which 
requires a different approach. We note that 


cos^ ^ Zg -I- sin^ A| = I, 


(A33) 


and use it in the following chain: 

||(A^5-Ab) 1^)11 

= 11(1-|Ab|) 1^)11 (A34) 

< 11(1+|Ab|)(I-|Ab|) 1^)11 (A35) 

= cot2(M)||(I-^i) 1^)11 . (A36) 

= cot2(^)||(I + Za^b)(I - ^a^b) 1^)11 (A37) 

< cot^(^)(l + (cos/r)“^)(5i = ^4 (A38) 


The first equality uses unitarity and the property 
AglAnl = Ab. The first inequality uses the operator 
inequality I + |Ab| ^ I- The se cond equality uses (A33|. 
The last inequality uses (A33) again to put a bound on 

PbIIoo- 

We now t urn t o the self-testing statement. In the isom¬ 
etry output (A4), because Z'a and have near-identical 
action over the state by bound (A32), the dichotomicity 
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of the regularized operators makes the two middle terms 
approximately vanish. In the first and last terms, for the 
same reason, the projectors (I ± and (I ± 

are nearly identical, and idempotence can be used. In the 
first term, the error introduced by this approximation is 
bounded as follows: 


I + Z'p^ I + Zq I + Z'^ 




{I + Z'j,){{I + Z{,)-{I + Z'^)) 


I^b-^aII 


<<5i. 


(A39) 

(A40) 


We used the fact that (I + = 2(1 + in the first 

equality, and the operator bound P+^aIIoo < 2 in the 
first inequality. By the same reasoning, and using the 
fact that and Ag are unitary to discard them from 
the norm, the fourth term leads to the same bound: 




l-Z'^l-Z^ I - ^A 


< ■ (A4I) 


The two middle terms are similarly bounded using the 
orthogonality of complementary projectors: 


I T I ± 
2 2 


-0 


<<5i. 


{ItZ'j,){{I±Z{,)-{I±Z'^)) 


(A42) 

(A43) 


Putting these bounds together, we deduce that replac¬ 
ing the isometry output with the state 

P) |00) + AiA(,^P^ p) 111) (A44) 


yields error terms bounded by 45i. We note that this 
state is only approximately normalized. 

We now take this approximation further to show that, 
as in the ideal case, the physical registers in the two 
terms are (approximately) proportional to the same 
state, which we call |junk) and define as 


ljunk) = /3 


_i I -I- Z'p^ 

2 cos Q 


(A45) 


where /3 > 0 is such that |||junk)|| = 1. W e will later 
show that /3 ~ 1. The first term in (A44) is therefore 


/I cos 0 ljunk) 100). We show proportionality for the sec¬ 
ond term: 


A' A'“"^a 


< 


— /3sind ljunk) 


X'kX'^ 


l-Z', 


J- 


— tan0- 


^A 


-^B ^ ^^^ -A^tand 


2 

I + ^B 


-f 


tan0 


^B-^A 


(A46) 


(A47) 


The first equality uses the definition of |junk). The in¬ 
equality uses the triangle inequality and the unitarity of 
A^ on the first term. The second term can be bounded 
by (A32). We can see that the first term is similar to 
bound (A20), from which we now derive an equivalent 


for the regularized operators: 

(cos6IA(5(I - Z'p^ - sin0A^(I -h Z'^)) p) 

^cos0Ab(I — ^a) — Aa sin 0(1 

+ ||cos 0(I-Z),)(A(5 -Ab)P)|| 

+ \\si^ex'^{z^-z'^) ” 


< 


Zb)) H) 


<82 + 284 cos 0 -I- i5i sin 0 = (Jn ■ 


(A48) 

(A49) 


The first inequality uses two triangle inequalities, the 
identities Z'a = Za and A^ = Aa and the commuta¬ 
tion of these oper ators with Bob’s. 

All in all, (A47) is thus bounded by tan 98\ + 82 ! cos 0, 
and we end up with the following bound showing an ap¬ 
proximate tensor product structure on the output state: 


$(|'!/j)) —/? ljunk) 1^/;) < (4-I-tan0)(5i 




2 cos 0 


= 5. 
(A50) 

The last error term to recover the state self-testing state¬ 
ment (A2) is to show that /3 ~ I. To do so, we bound 


1 — /3 from above and below. First, using the property 
that the isometry preserves the norm, we find 

1 = ||$(|^))|| (A51) 

< II^’(1^)) -/? ljunk) |V^)|| -h/3|| ljunk) |V’)|| (A52) 

<8 + 13. (A53) 

Similarly, we have 

11/3 ljunk) It/;) II (A54) 

< 11/3 ljunk) \iP) - $(|VI))|| -f ||«>(|^/i))|| (A55) 

<8 + 1. (A56) 

Put toge ther, /3 is constrained as \1 — (3\ < 8 . Using 
this and (A50), we finally complete the proof for the self¬ 


testing statement on the state 


^(1^)) - ljunk) 


< 28. 


(A57) 


h. Self-testing the measurements 


The proof for the self-testing statements involving op¬ 
erator actions in (A2) builds on the previous section. As 


a first step, we give approximations on the output of the 
isometry when acting on the input with the regularized 
operators, and then further approximate by substituting 
them with the nonregularized operators. 

This is where we need the approximate anticommuta¬ 
tion of and Z'a when acting on the state, given by 
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inequality (A22). We will also need to prove a similar 
anticommutation bound for and As we noted 
earlier, the nonregularized Xb and Zb anticommute by 
definition, but the regularization procedure breaks this. 
As a result, we can only prove that {Xb,^b}|V') — Oj 
which will require a number of successive approximations. 

The use of anticommutation is made clear when look¬ 
ing at the action of the isometry (A4) on X'j^\ij}) or 
Ag \tjj). For example, (the reasoning for the other op¬ 
erator is the same): 


1^)) 


= -[{! +z'B){i-x'^z'M\^)\m} 

+ Aa(I -I- Zb){ 1 + X'p^Z'p^X'p^ IVj) |00) 

^ x'^{i- z'^){\- x'^z'Mm^^) 

+ x'M^ - z'b){i + m |01)]. 

(A58) 


The approximate anticommutation of and means 
that X'p^Z'pX'p^ |'(/j) ~ —Z'p^ IV^), and this output only dif¬ 
fers from $(!'!/;)) by the action of Ox on the first ancilla 
qubit, and a small error term. On the other hand, it is 
easily seen that the isometry applied on Z'^ \i,) or Z'^ 1^) 
yields exactly the output d)(|^/!)) with the cr^ Pauli oper¬ 
ator acting on the corresponding ancilla. 

We now prove the anticommutation bound for Bob’s 
operators. Compared to the other bounds we have de¬ 
rived so far, this is not as immediate as it might seem; 
indeed, while we know by ( |A30 ) and (A381 that the non¬ 
regularized and regularized A and Z operators on Bob’s 
side are approximately interchangeable, this is only the 
case when they act on the physical state j'i/;). Thus, 
we can not deduce from this alone that {Ag, Zi,} M ^ 
{Ag,Zg}|V') = 0 because not all unitaries in the anti¬ 
commutator act directly on the state. Instead, what we 
do is to approximate the rightmost operator in each term 
of the anticommutator by its unnormalized counterpart 
by using (A30l and (A38), and translate the action of 


Zb and Ab into one that commutes with the leftmost 


operators using respectively (A19) and (A21|. Then, the 


leftmost operators can in turn be approximated because 
they now act on the state directly. 

The steps outlined above are carried out as follows: 

II{^b,^b} 1^)11 


<2S^ + \\iZA + Zi,)X{,M\ 

(A59) 

< 26i + 2^4 -t 2(55 


+ + ^b)^a(I - cZb) 1^)11 

(A60) 

= 26i + 25i + 2(55 


+ s-i(I - cZb)(^A + ^b)^A 1^)11 

(A6I) 

< 26i + 25i + 2(55 


+ l^^'/cos ^11 ZaAa+AaZ(5 1^)11 

(A62) 


< 25i + 2^4 + 2(55 + 

^<5?. 


1 + c/cos /i 




(A63) 

(A64) 


In the first i nequa lity, we used ( A32 ). The second in¬ 
equality uses (A38l followed by (A21), which lets us com¬ 


mute this approximation of Ag to the left of the operator 
product in the first equality. Next, in the third inequal¬ 
ity, the operator bound ||^b|Ioo — (cos/r)“^ that we used 


to derive (A38) is used again. Finally, we use (A30) fol¬ 


lowed by Alice’s anticommutation bound (A22) to reach 
the last inequality. 

Thus, combining the anticommutation bounds with 
the regularization approximations for Bob’s operators 
( A30|) and (|A38|), we find 


$(Za|^))-(7^<I>(|^))|| =0, 

(A65) 

$(AA|^))-aA$(|^))|| <25^ 

(A66) 

$(Zb|^))-(7MI^))II <<5i, 

(A67) 

$(Ab|^))-(t«4>(|^))|| <(54 + 2,5^ 

(A68) 


In the last two bounds, we have also used the fact that 4) 
preserves the norms and is linear, such that for example 

mzBm-^{z'^\m = \\{ZB-z'^)m. 

Our goal is to reach bounds from the joint action of 
the observables Ax and By. So far, we can only compute 
bounds for the action of one party at a time. Indee d, to 
bound W^iBy \'ij;)) — By^{\’tp))\\ we can use definition (A1) 
for By, the triangle inequality, and the bounds above for 
the action of Zb and Ab. 

We now show that with joint action of both parties, 
Alice’s operator is easily dealt with. First, consider 
4>(AoBy As with (A65), this is exactly the same 
as AQ^{By I'i/j)) because Aq = Z'p^ and Aq = a^. Next 
is ^{AiBy I'ljj)), for which it is easily seen that the rea¬ 
soning that we used to reach (A66) is unchanged. Indeed 
Ai = A^, and this state is identical to (A58) with \ijj) re¬ 
placed with By\tp). We use commutation between Alice 
and Bob to move By to the left of (I ± X'pZ'pX'p^ and, 


because \\B. 


y iloo 


= 1, we find 


mA,By\^))-AMBym\<2Si 


(A69) 


We have thus showed that (A65) and (A66) are un¬ 
changed if we replace \tjj) with By \tjj) in both terms. 

In addition to the above, we will use one last ap¬ 
proximation to replace the isometry output 4>(|'0)) with 
ljunk) |'0) , at the cost of the additional error term of 
26 from (A57). We finally find the following bounds 
Cxyie) in the self-testing statement (A2|, where we de- 
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fine €xy 


— ^'xy + 2^: 



= ®0,-l 

= 0, 

(A70) 



(ATI) 

p' — A 

^- 1.0 — ®-i,i 

— A —A 

— ^ 0,0 — ^ 0,1 


= cos 

fj.+ {64 + 2(5®) sin/r , 

(A72) 

0 

II 

II 

: 26^ + e'_i Q . 

(A73) 


Appendix B: Additional SOS for CHSH 


By guessing values for the parameters in ( fid]), o ther 
nontrivial SOS decompositions than ( [T7| ) and (|18[) can 
be found. 

We report here a combination that leads to a SOS ma¬ 
trix of rank 4 which is extremal in the set represented in 
Figure i.e., it cannot be decomposed as a convex com¬ 
bination of SOS matrices of the same form as (141. With 
the values q = 1/4, = 5/8, A = 1/4, the square root 

of M is rather well-behaved, and leads to the following 
SOS decomposition: 


^0 = 


1 

8^2 


2{Za - Z^f + 5(Aa - Xb)2 
+ 2(ZaAb+Xa^b)" 


+ (3 - 2ZaZb - XaX^Y 


(Bl) 


An additional symmetry of the CHSH inequality can 
be exploited to reach an additional SOS decomposition: 
swapping X and Z (i.e., swapping the observables Aq 
and Ai and changing Bi to —Bi) leaves Iq invariant, 
and therefore 


Io = 


1 

8^2 


2(Aa - Xsf + 5 (Za - Zsf 
+ 2 (ZaXb + XaZb)^ 


+ (3 - 2AaAb - ^a^b)" , 


(B2) 


for which the parameter values are q = 1/4, /i = 1/4, 
A = 5/8. 

It is not clear whether either of these decompositions 
yields an easy generalization to the whole la family. 


Appendix C: Shortcomings in previous results 


1. CHSH self-test 


In their robust self-test proof for the maximal CHSH 
violation, McKague et al. [12] introduce the isometry 
used in [7| and in the present article, with the same reg¬ 
ularization construct we used on operators (A3) to build 


the unitaries that make up the isometry. The proof tech¬ 
nique is the same as in this article: they identify from 


the self-testing criterion a series of constraints on the ac¬ 
tion of the observables in the system, that are then com¬ 
bined to form the self-testing statement. However, in 
their self-test of the observables, the authors rely on the 
anticommutation relation {Ag, on Bob’s regu¬ 

larized operators. As we noted, although Xb and Zb (in 
our notation) anticommute by definition, this property 
is lost to regularization. Indeed, the nonzero eigenval¬ 
ues of anticommuting operators come in pairs of oppo¬ 
site sign, which means that nonsingular anticommuting 
operators do not exist in odd-dimensional Hilbert spaces. 
In fact, if Bq and Hi do not each have as many -|-I as 
— 1 eigenvalues, their eigenspaces must share a nontriv¬ 
ial intersection, which will be an eigenspace for both Ag 
and (with nonzero eigenvalue because of regulariza¬ 
tion) where these two operators can therefore not anti¬ 
commute. Hence, in their proof in Appendix B of [12] . 
while Bq ± B'l anticommute, it is not true anymore when 
their zero eigenvalues are replaced with I. 


However, this oversight only affects the final result by 
error terms of the same robustness order as they claim. 
Indeed, as they show in the proof for their second self¬ 
testing criterion based on Mayers and Yao’s work, the 
anticommutation of Ag and ^B in front of the physical 
state |'0') can still be given a robustness bound. This 
follows from the O(^) bound on the anticommutation 
of Alice’s A^ and Za operators and the 0(e^/'^) bounds 
on the replacement of Bob’s A' and Z' operators by Al¬ 
ice’s, which can be combined to transform Alice’s anti¬ 
commutation bound to one for Bob with 0(e^/^) order. 
Although this is worse than the £2 G 0(^/e) that they 
use in their Theorem I, the final order in the self-testing 
bounds is unchanged by this correction because they al¬ 
ready contain 0(e^/^) terms from Alice’s anticommuta¬ 
tion bound. 


2. Partially entangled state self-test 


In their proof for the self-test statement (A2) on the 
state (i.e. x = y = —1), Yang and Navascues introduce a 
SOS decomposition for different from ours. |7] They 
write the SOS polynomials Pi as products q.^ • V with V 
defined as in (|^ , and decomposed in terms of five nine¬ 
dimensional vectors different from the ones we defined 
in We reproduce here the qi they list in their article, 
with a change in 52 which originally contained a typo that 
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was communicated to us m- 

qi = 




1/2571 + 52-9-778 

q2 = -^——-— (i-i + era - erg) 

10s 

27 - 25c73^ 3 , 

qa =--ri + T7r(r5 - r4), 


3072 


10 ' 


35 , , 5 c 73^ 

q4 = 3^(r, + r,)-^j^r,, 

7 4972 + 9800 c 773 — c + 7 
q"; = ---ri , 

45 ^20 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 


with c = cos(40); 7 = y (75 + 25c)76 — 2c — 72; and 
w = 18125 cos( 86 ») - 72500 cos(46i) - 108706, Their SOS 
is then defined as 




(C 6 ) 


This decomposition is problematic. Indeed, the vec¬ 
tor qs is real only in the interval in 6 (or equivalently 
in a) where the expression inside the square root in 
(C5) is positive. This corresponds approximately to 


] 0 , 7 r/ 4 ] = ]0,0.78540]. Outsi de o f that interval, a change 
of sign of the fifth term in (061 is actually required to 


recover the left-hand side. However, this change means 
that the decomposition of I a, while valid, is not a sum of 
squares anymore, and therefore it does not have the prop¬ 
erties required for the self-test proof. The value 9 = 7 r /8 
is within the validity interval; we represent the SOS in 
Figure]^ as a red dot. 

Another issue with this decomposition is that the 
vectors have a linear dependency that went unnoticed in 
the original article. This can be seen in the fact that 
all five qi only depend on r 4 and through their dif¬ 
ference rs — r 4 , which means that the qi only span a 
four-dimensional subspace. Hence, this is insufficient to 
certify the five operator identities in Si+ab- Notab ly, 
forming the operator on the left hand side of (A16) in 


9 G [0.07574,0.73014], which is not the full interval 


the present article from their polynomials qi • V requires 
a decomposition of rs + r 4 in terms of the qi vectors as 
noted in the supplemental information to their article, 
which is not possible. 

The linear dependency between the qi is also visible on 
Figure]^ Indeed, the points on the boundary of this set 
correspond to singular SOS matrices M, as the nonsingu¬ 
lar matrices, being strictly positive definite, do not satu¬ 
rate the inequalities that define the boundary. Therefore, 
for this SOS decomposition, M is singular (i.e. of rank 
at most 4 here) and the qi can therefore not be linearly 
independent. 
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TABLE I. Vertex SOS decompositions for CHSH. 

Cl \ — 0 , fi — 0 , q = 0 

/O 0 0 0 0\ 

, [ 0 1 0 0 0 

M = —^ 0 0 1 0 0 , (C7) 

2^2 0 0 0 0 0 
Vo 0 0 0 0/ 

^ 2^ ^ ^a^b)" + (-1 + VaXb)"] ; (C8) 

C2 A = 1, /i = 0, g = 0 

/I 0 0 0 0\ 

, [00000 

M = —- 0 0 1 0 0 , (C 9 ) 

2^ 00000 
Vo 0 0 0 0/ 

To = ^ [Rl +Rl] = -^ {{-Za + Z^f + (-1 + XaXb)"] ; (CIO) 

C3 A = 0, /i = 1, (j = 0 

/O 0 0 0 0\ 

, [01000 

M = 0 0 0 0 0 , (Cll) 

2v2 0 0 0 1 0 

Vo 0 0 0 0/ 

To = 2;^ [Rl + ^4] = ;^ [(-I + Ta^b)" + (-Va + Xb)"] ; (C12) 

C4 A = 1, /i = 1, g = 0, 

/I 0 0 0 0\ 

, [00000 

M = —- 0 0 0 0 0 , (C13) 

2^2 0 0 0 1 0 

Vo 0 0 0 0/ 

To = ^ [Rl +Rl] = ^ [{-Za + Zsf + {-Xa + Xb)"] ; (C 14) 

Cs A = 0, M = 0, g = 1/2 

/O 0 0 0 0 \ 

, [01/2 1/2 0 0 

M=— 0 1/2 1/2 0 0 , (CIS) 

2v2 0 0 0 0 0 

Vo 0 0 0 1/2/ 

To = 4;^ [{R2 + Raf + ^5] = 4:^ [{-2V2 I + To)" + 2 (ZaXb + XaXb)"] + ’ (^16) 

where S' = Ao(-Bo — Bi) + Ai(i3o + Bi)- 



